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Abstract 

In the presence of quantum measurements with direct photon detection the evolution 
of open quantum systems is usually described by stochastic master equations with jumps. 
Hcuristically, diffusion models can be obtained from these equations as approximation. A 
condition for a general diffusion approximation for jump master equations is presented. 
This approximation is rigorously proved by using techniques for Markov processes which 
are based upon the convergence of Markov generators and martingale problems. This 
result is illustrated by rigorously obtaining the diffusion approximation for homodyne 
and heterodyne detection. 
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Introduction 



In open quantum systems [H HU HH [El [231 H9], especially in the theory of measurement, an 
active line of research concentrates on the study of stochastic master equations and stochastic 
Schrddinger equations These 
equations are stochastic differential equations (SDEs) describing the evolution of a small 
system undergoing continuous measurement. Specifically, stochastic Schrddinger equation 
describe the evolution of vectors in Hilbert space {wave function of the small system) [T3"l 
fTH| [3 El]- Moreover, the notion of stochastic Schrddinger equation is closely related to the 
concept of unravelling in quantum mechanics which is highly used to derive Monte Carlo wave 
function simulation methods (see Chapter 7 in [3] and Chapters 17-19 in |15j for examples, 
applications and for other references). The term stochastic master equation refers to the 
evolution of density matrices (stated of the system) [3 HI El [351 EH EH CEB] (see [3 El [351 [36] 
for the equivalence between these two types of equations). 

The framework sketched above has a wide range of applications in quantum optics and 
quantum communications (SJISIiaEltTnillllliaEaiia A 
particularly important application is the modeling of photon detection. A typical situation 
is the one of a small system, e.g., a two level atom, stimulated by a laser. The measurement 
scheme is the continuous detection of the fluorescence photons emitted by the two level system. 
In this situation, the evolution of the state of the system is given by the following stochastic 
master equation 

d Pt = L( Pt )dt + ( T ^Cp t C*] ~ P ) ^ " rTc [ C Pt C *] dt )> (!) 
where (JV$) is a counting process with stochastic intensity JqTic[C p s C*]ds. In this equation 



1 A state or a density matrix is a positive trace class operator of trace one. 



the operator L is called Lindblad operator and takes the form 

L(p) = -i[H, p] - ±{C*C, p] + C P C\ (2) 

where C is an arbitrary operator. The Lindblad operator is actually the generator of the 
reduced dynamic of the qubit if there were no measurements. The occurrence of a jump in 
equation (pQ) corresponds to the emission of a photon by the qubit. Another type of equation 
describes diffusive evolutions, that is, 

d Pt = L( Pt )dt + (Cp t + Pt C* - Tr[p t (C + C*)])dW u (3) 

where (Wt) is a standard Brownian motion. Such an equation is related to heterodyne or 
homodyne detection (the first one allows to study the spectrum of the light and the second 
one the squeezing [21 El QSl El CDS HQl H2] ) • 

Heuristically, equation (|3|) can be obtained from the model ([T]) by the following approach 
[21 H31 H51 EH 1251 H3 H2]- In equation CP), we replace the operator C by C £ = C + I/e, and 
we define the process (Wf) by dWf = edN t — dt/e. In this way, formally by using usual Ito 
stochastic calculus we get dW^dWf = e 2 dN t = edWf + dt and by taking the limit e goes to 
zero, we get lim e _>oWt £ = Wt, where Wt is a Brownian motion. Next, by assuming C = C* 
and by taking the formal limit, e goes to zero, in 

dp t = L e [p t )dt + \ ( £[c%*\ ~ Pi ) ^ -^[C'PtGW' W 

we recover the equation §3§ . Such a result is called diffusion approximation (see |16[ [T71 [BUI [55] 
for similar result in biology, classical mechanic or mathematical finance). Let us stress that 
the equations (pQ) and ([3]) are particular cases of stochastic master equations. More general 
situations are described by jump-diffusion SDEs [7[ [Ml [37]. 



3 



The main aim of this article is to prove rigorously the diffusion approximation. We 
address this question in the more general models of jump-diffusion SDEs. We present a 
general condition for the coefficients defining the SDEs which ensures diffusion approximation 
results. We show that this condition appears naturally in concrete physical applications. 

The article is organized as follows. In Section 1, we investigate the diffusion approxima- 
tion problem for the jump-diffusion stochastic master equations. The approach is based on 
convergence of Markov generators and problems of martingale. Then, we obtain a sufficient 
condition where the diffusion approximation is valid. In Section 2, we show in concrete exam- 
ples, that such situations appear naturally. In particular, we justify rigorously the diffusive 
models of homodyne and heterodyne detection in quantum optics. 

1 Diffusion Approximation of Jump-Diffusion Stochastic Mas- 
ter Equations 

1.1 Framework and Main Result 

We consider a Hilbert space 7i = C , which represents the small system, which undergoes 
indirect quantum measurement. In order to describe the random evolution of the small sys- 
tem, we consider a probability space (Q,J-,J-t,P), which supports a n dimensional Brownian 
motion (W t = Wi(t), . . . , W n (t)) and k independent Poisson point processes Ni(),i = 1, . . . , k 
on IR 2 , independent from the Brownian motion. Let us remind that, for all Borel sets A of M 2 , 
the quantity Ni(A) corresponds to the number of points of the Poisson point process Ni into 
the set A. Moreover, for each i, the application A i — ► Ni(A) defined on the Borel sets defines 
a measure. More precisely, this defines a Random measure |27j . Moreover, for all Borel sets 
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A, we have E[iVj(.A)] = X(A), where A corresponds to the Lebesgue measure on M? (here, 

the Lebesgue measure is called intensity measure of the Poisson point processes [27]). Now, 

by denoting the differential elements as Ni(dx,ds), we can consider the following stochastic 

differential equation 

ft n rt 

Pt = ^ L(p s .)ds + J2j o (c iPs _ + Ps _C*-Tr[p s .(C + C*)]yW i (s) 

+ ^2j Q J R I Tr [d'p D*} ~ Ps ~) lo < x < T ^ D 'Ps-D*]( N i( dx ' ds ) -Tr[Dip s .Di]dxds), 

(5) 

where L is a Lindblad operator, defined by 

1 n i ^ n k 

L(p) = -i[H, p]-- Y,{CtC uP ] - - J2{DtDi,p} + C lP Ct + D iP D*. (6) 

i=l i=l i=l i=l 

In the equation ©, the operators Cj, i = 1, . . . , n and Di, i = 1, . . . , k are N x N matrices. 
The operator H is a self adjoint operator called free Hamiltonian of H. The equation (0) is 
a generic stochastic master equation (see [7J[31] for different versions). Let us stress that, by 
using a Poisson point process on M 2 , we can define rigorously the counting processes: 

Ni(t)= / / l 0<x< Tr[D iPa .D*]Ni(dx,ds), i = l,...,k, 
Jo Jr 

which are clearly counting processes with stochastic intensities Tr{Dip s _D*]ds. The solution 
of ([5]) is called a quantum trajectory and takes values in the set of states of H; this describes 
the random evolution of Ti. in presence of continuous indirect measurement (see Refs. [7J [Ml 
[35J1361E7J for justifications on existence and uniqueness of solutions). In the rest of the paper 
we denote S the set of states, this is a compact subset of B(TL). 

Our aim consists in generalizing the heuristic rule described in the Introduction and to 
prove rigorously diffusion approximations in such models. The mathematical framework is 
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the one of the convergence in distribution for stochastic processes. To this end, let T > 
be fixed, we introduce the space T)[Q, T[, which corresponds to the set of cadlag processes 
endowed with the Skorohod Topology (topology of convergence in distribution [221 l38j 125])- 
The main result is expressed in the following theorem. 

Theorem 1 Let J be a subset of I = {1, . . . , k}, let define for all j £ J, D £ - = Aj+I/e, where 
Aj,j S J are N x N matrices. Let Ni(t) , i = 1, . . . , k be k independent Possion point processes 
on R 2 , let (Wt = (Wj{t))j & j) be a car d(J)- dimensional Brownian motion independent of the 
Poisson point processes. 

Let assume that the condition 

E " 4) = (7) 

jeJ 

is satisfied. Therefore, the solution of 



p\= f (l(pI)+ ]T (-D iP iD* + Tr[D iP iD*}pi) 
Jo V i£l\J 

+ ^ [-D*p{D*y + Tr[D E jP E s XDj)*]p e s _) )ds 




converges, when e goes to zero, in the space T>[0, T[, for all T, to the solution of the stochastic 
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differential equation 




It is then straightforward that the above Theorem is a generalization of the heuristic 
approach presented in the Introduction. The next subsection is devoted to the proof of this 
result. 

1.2 Proof of Theorem 1 

The proof of Theorem 1 relies on two steps. First we prove that, under the condition (J7|), 
the family of processes (pf) owns the tightness property which corresponds to the relative 
compactness criterion in the Skorohod Topology. Secondly, we show that the family of Markov 
generators (A e ) associated with the processes (p £ t ) converges to the Markov generator of (pt). 
Next, by combining the tightness property with this result, we get the final convergence. The 
tightness property follows from the next proposition. 

Proposition 1 Let T > be fixed. Let assume that the condition ([?[) is satisfied. Let (pi) be 
the solution of the equation §3$. For all M > 0, there exists some constant Z such that for 
all e < M 

v[M 2 -M 2 M-p £ tl \\ 2 ] <z(t 2 - tl ) 2 , (io) 
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for all t\ < t < t 2 < T. 

Therefore, the family of processes (pt)o<t<T,£ > is tight for the Skorohod topology on 
V[0,T[. 



Proof: The fact that the property (jlQI) implies the tightness for the Skorohod topology is a 
classical result (see Theorem 15.6 in |12] and Ref [29J for further explanations). Here, we just 
show that the inequality (|1U|) is true. Before attacking the estimation, we need to notice the 
following two facts. 

First, since the condition ([7]) is satisfied, an easy computation gives 

L e (p) = L(p), (11) 

for all states p £ S. As S is compact, the function L e is bounded by a constant K on S 
independently of e. 

The second fact concerns the estimation of the terms D e -p{D^Y /Ti[D E -p{D E -)*\, j 6 J. For 
all j G J, with the definition of D £ , we get 



p + e[A 3 p + pA) - Tr[p(Aj + A*)]p) + o(e). (12) 



D £ 3 p{Dlf 
Ti[DCp(Dt)*} 

It is important to notice that the o(e) are uniform on S (since S is compact). 

Now, we are in the position to prove the inequality (|10|) . To this end, we introduce the 
natural filtration (Qf) of (pf), that is, = a{p^,u < t}. In this way, we first estimate the 



S 



term A(t,t 2 ,e) = B[\\pf 2 - p% \\ 2 \Gf}. We have 
A(t,t 2 ,e) 



< E 



i£l\J 

+ E E 



\\m.)\\ds 

t2 



Apia* 

Tr[ApiA 
D*pi(D 



Pi- ) 1 0<x<Tr[D I pl_D*]^(d2;,(is) 



Ps- !o<:r<Tr[iXpl_(Df)*]Nj(<ix,cZs) 



5? 



(13) 



Let (i), (ii), (Hi) denote the three above terms of the sum. 

Since L is bounded on 5, for the first term, we have (i) < K(t 2 - t) 2 < KT(t 2 - 1) almost 
surely. 

By remarking that card{I \ J} < k, for the term (u), we have 



(ii) < k sup E 
ie/\J 

< k sup E 

iei\J 



TY[Apf D*} ~ Ps -J i °< Tr [ApI.O*] 1 0< :C <Tr[D,p|_D* 



2 




dxds 





Apia 



Pl- ) 1 0<Tr[D I p|.D*; 



Tr[D iP iDt]ds 



0! 



(14) 



Tr[ApI.A 

The passage from the expression (ii) in (I13h to the first inequality in (|14|) relies on two 
facts. First, we have computed the expectation with respect to the Poisson point processes 
by using the usual property of stochastic integration theory with respect to random measures 
(such a property is called Ito isometry property). In other words we have used that the 
intensity measure of each Poison point process Ni is the Lebesgue measure (we refer to |27] 
for a complete introduction of stochastic integration theory with respect to random measure). 
Secondly, we have introduced the term lo<Tr[Dip| D*] by remarking that lo<a;<Tr[_Dip| D*] = 
lo<Tr[Dip| D*]lo<x<Tr[Dip% D*] (this allows to keep the fact that we have the property < 
Tr[Dip e s _D*] in the second inequality). Thus, with respect to the underlying probability, it is 
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important to notice that we implicitly consider that the quantity Dip e s _D* /Tr[DiPg_D*] = 
if Tr[Dip £ s D*] = 0, otherwise this quantity defines a state. Now, since the set of states S 
is compact, the term inside the L2 norm is almost surely bounded, then the term inside 
the integral f^ 2 is bounded. As a consequence, there exists a constant K\ such that (ii) < 
K\ (£2 — t) almost surely. 

For the term (Hi), since card{ J} < k, by using the estimation (I12p . we have 



(in) < ZcsupE 
j'eJ 



< fesupE 



D*pi(D 



*2 



pi 



P%- 1 0<s<Tr[DfpI„(Df)* 



dxds 



Qt 



Qt 



< fcsupE 



" + -TrlpKAj + + o( £ )" 2 



i + ~Tr[pt(^ + 4)] + TtLV^A ds 



6? 



(15) 



Now, it is straightforward to notice that the factor e, in the L2 norm, compensates the 
terms 1/e 2 and 1/e. Hence, for M being fixed, there exists a constant K2 such that for all 
e < M, the term inside the expectation is bounded (recall that S is compact). Then, we have 
(Hi) < .^2(^2 — t) almost surely. As a consequence, for all M > 0, there exists a constant 
S = KT + K\ + K2 independent of t, £2 and e, such that for all e < M, we have almost surely 

A(t,t 2 ,e) < S(t 2 -t). 

We shall show that this implies the expected result. By conditioning with Q\, we have 



^[M 2 -Piw 2 \\Pt- P £ t : m 



.£ I|21 



E[A(t,t 2 ,e)M-pl 



< S(t2-t)E[A(t u t,£)] 



< S 2 (t 2 -t)(t-h) 

< ^(t2-h) 2 



(16) 
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and the result is proved with Z = S 2 /2. 



□ 



Now, we address the second step which concerns the convergence of Markov generators. 
Moreover, this convergence gives a reverse result concerning the condition ([7]). Before ex- 
pressing the result, we introduce the notation C 2 for denoting the set of C 2 functions defined 
on BiTL) with values in R and with compact support. The terms, D p f(.) and D^f(., .) will 
denote the first and second differential. Besides, we define the functions hj,j G J on S by 

hj{p) = A jP + pA* - Tr[p(Aj + A*)]p, 

for all states p. 

Proposition 2 Let A e be the Markov generator of the process (pf) defined by 



AJ(fi) = D„f(L e (p)) 



E 

ie/y 

E 



f 



f 



Tr[Dlp(DI) 



Hp) - D„f 
- Hp) - D„f 



p l pp; 

Tr[D lP D 
( D*p(D 



3 > 



Tr[D iP Dt] 

Tr[D £ jP (D 



f)1 



(17) 



for all / G and for all states p. 

Let A be the Markov generator of the process (pt) defined by 



Af(p) = D p f(L(p)) 



Tr[D iP D* 



+ E / 

iei\J L 

+\Y, D2 P fi h MMp)), 



Hp) - D p f 



DjpDj 
Tr [DipD* 



Tr[D iP D* 



(18) 



for all f £ C 2 and for all states p. 
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Then, let f G C%, we have the convergence 



Urn l sup \A £ f(p)-Af(p)\ = 0, (19) 
pes 



if and only if the condition Q: Ylj^j y^j ~ A*j = is satisfied. 
As a consequence, we have 



limE 

e->0 



/t+s \ »" 

' i=l 



0, (20) 



/or a// m > 0, for all < t\ < ti < . . . < t m < t < t + s, for all functions {9i)i=x t ,.. m and for 
all f in Cg. 

Remark: The fact that the operators A e and are the Markov generators of (pf ) and (pt) 

follows from Ito stochastic calculus (see [22, 27J concerning the definition of Markov generators 

and [37J for explicit computations). Moreover, we have the following important property. 

Let po be a state and let / € C%, the process defined by 

ft 



f(pt) ~ f(Po) ~ [ Af(pt)ds, (21) 
J o 



for all t > 0, is a martingale with respect to the natural filtration (Qf) of (pf), where £/f = 
o~{pf l ,u < t}. This property is related to the notion of Martingale problem. Let us make precise 
this notion. In probability theory, the couple (A £ ,po) defines what is called a Martingale 
problem. Solving this martingale problem consists in finding a Markov process (Xt) such that 
Xq = po and such that, for all / G C^, the process (f(X t ) — f(Xo) — J* * A e f(X s )ds) is a 
martingale with respect to the natural filtration generated by (Xt) (see |27} [22"1 [29] for more 
general definitions of martingale problem). In this sense the process (pf) is a solution of the 
martingale problem (A e ,Po) (the same holds for (p t ) with (A, Po))- 

Moreover, as (pf ) is the unique solution of the equation ([8]) and as the property (f2Tj) is 
satisfied, the process (pf) is the unique solution, in law, of the martingale problem (A e ,po) 
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(see Refs. \22\ [29], [28] |27] for general considerations on martingale problems and [37] for 
results concerning stochastic master equations). Note that the process (pt) is also the unique 
solution, in law, of the martingale problem associated with (A, po). 

Since we shall prove a result of convergence in distribution, the result of uniqueness, in 
law, is crucial. Indeed, if the limit process of (pf ) satisfies the martingale property (f2T1) for 
(A,po), this limit process will be equal to (pt) in law. Let us prove, now, the proposition. 
Proof: Let j G J, the estimation (fT2"j) stated in the proof of Proposition [U allows to apply 
the Taylor formula. This gives 



/ 



Tr[D £ p(D £ y 



Hp) - D p f 



Tr[Dtp(Dz)*] 
1 



Tr[D £ jP (D 



|)1 



D 2 J{h 3 (p),h 3 {p))+o(e), 



(22) 



where o(e) is uniform in p. As a consequence, it is easy to see that 



lim sup 



pes 



E 

.jeJ 



Tr[D £ p(D 



Hp) - D p f 



DjP(DjY 
Tr[D £ iP (D 



Ti[D £ p(D 



0. 



(23) 



This first approximation shows how the diffusive part, in terms of Markov generators, appears 
in the limit. Let us treat now the term D p f(L £ ). The only contribution involves the terms 
where the operators D £ - appear. For all j G J, we have 
1 



D p f [ --{(D £ YD £ ,p} + D £ p(D £ y 



D p f (-±{AjAj,p} + A jP A^j + L Dpf (( Aj - A*)p + p{A j - A*)) . 



(24) 



As a consequence, the two expressions (|23|) and (|24j) imply the equivalence 



lim sup \A £ f(p)-Af(p)\ = 0,V/ G Cl & Yl { A i- A i)p+p{ A i- A i) = °^P G S - ( 25 ) 



pes 
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Now, if we have YljeJ ( (A? ~~ A^)p + p(Aj ~~ ^J) ) = 0, for all states, this equality holds for all 



positive matrices (by multiplying by the trace). Since any Hermitian matrix can be written 
as linear combination of positive matrices (this comes from the diagonalization) , this equality 
holds for all Hermitian matrices. Now, since we have X = (X + X*) /2 + i(—iX + iX*)/2, for 
all matrices X and since the matrices X + X* and —iX + iX* are Hermitian, the equality 
holds for all matrices. Thus, we get the complete characterization ([7]) and the equivalence is 
proved (actually the uniform convergence is not necessary to imply the condition ((7J)). 

In order to finish the proof of the proposition it remains to prove the convergence (I20p . To 
this end, we insert the term f^ +s A £ f(p £ u )du in the expectation. Next, according to the remark 
before the proof, we use the fact that f(pf) — f(po) — f A e f{p £ u )du is a martingale with respect 
to the natural filtration (Qf) of (pf) (we have to notice that Y\^=i @i{Pt ) ^ s @t measurable) . 
Next, the uniform convergence result of Markov generators ensures the convergence (120(1 . □ 

Now, we are in the position to prove Theorem 1. Let us assume that the condition ([7j) is 
satisfied. Let (e n ) be a sequence converging to 0. Since (pi") is a tight sequence (Proposition 
1), we can extract a convergence subsequence. Let denote (pt) the limit. According to the 
convergence (f2"U|) . the process (pt) satisfies 



for all m > 0, for all < t\ < t2 < ■ ■ ■ < t m < t < t-\-s, for all functions (0i)j=i,..., m and for all 
/ in C%. This implies that (pt) is the solution of the problem of martingale associated with 
A. By uniqueness, in law, of the solution of the martingale problem, the process (pt) has the 
same distribution than (pt)- In this way, every convergent sequence of processes (p £ t n ), where 
(e n ) converges to 0, converges in distribution to (pt). As a conclusion, the family of processes 
(pf) converges in distribution to (pt). 




i=l 
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Figure 1: Experimental setup of homodyning. 



D 2 



H L {t) 




local oscillator 



2 Applications 

This section is devoted to present situations where the condition ([7]) is naturally satisfied 
and where the diffusion approximation can be performed. As a first direct application, we see 
that the heuristic approach presented in the Introduction is rigorously justified by Theorem 1. 
Indeed the condition C = C* is no more than the condition ([7]). In the two following subsec- 
tions, we investigate the diffusion approximation in the model of homodyne and heterodyne 
detection. 

2.1 Diffusion Approximation in Homodyne Setup 

The model of homodyne detection describes a two-level atom whose emitted light is detected 
by photodetection. The source, i.e, the two level atom, is driven by an external interaction. 
The light, emitted by the atom traverses a beam splitter and interferes with the light field of 
a local oscillator. Then, two detectors D\ and D 2 detect the light (see Section 6.4 of |18| or 
[31 [4"0l W\\ Fi2l [T4"l [T5] for more details on the experimental setup and optical considerations, 
see Figure [[]). Typically, we record two different types of jumps according of which detectors 
detects the light. 
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Our aim is not to derive the equation describing this situation, we refer to [21 El El [T3J [THJ 
EH EH E51 GS EH EH E2] for justifications. Let us give the stochastic master equations govern- 
ing the evolution of the two level atom. To this end, we consider the stochastic Schrodinger 
equation, that is, the stochastic differential equation which describes the evolution of the 
two- level atom in terms of wave functions (or pure states). Here, the corresponding equation 
is given by 



o 



iH = "00+ / -i[H+ —\\(T.ip s -\\ )^s-ds 



\ \\{a. + i0)M ~**) 1 o<-<llK-^)^ll 2iVl(dx ' ds) 
where N\ and N2 are two Poisson point processes on M 2 . 

Remark: In this equation, the constant 70 represents the spontaneous emission rate. As we 
can see in the equation (|26p . this parameter appears in the stochastic intensity. The quantity 



/3 represents the amplitude of the local oscillator. The operator H is the free Hamiltonian 



of the two-level atom and <r_ 



matrix. 



1 



is the usual notation for the corresponding Pauli 



The corresponding equation for (pt), in terms of density matrices is obtained by defining 
Pt = IV'tK'^l (t ne notation \ipt){^t \ corresponds to the Dirac notation for the one dimensional 
projector on Cipt] one can find also the term pure statqj to define such kind of state). The 
equation for (p t ) can be derived using the Ito rules (see O El Ell ES] for computations). We 



2 A pure state is a state which is a one-dimensional projector. 
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get then the following stochastic master equation 

Pt = Po + f Q (-i[H,Psl-\{(<T-+iP)*((T- + if3),p s -}-^{(v--iP)*(a.-iP),p s ^ 

+ Tr[((j. + i(3)p s _(a. + i/3)*]p 8 . + Tr[(a. - i(3)p s _(a. - i/3)> s .)ds 

f* f ( (o-. + i/3)p s (a. + i/3)* \ 
+ Jo L \Tr[{o-+iP)p s _{o_ + l (3y}~ Ps -) 1 o<^<?T r [(.. + ^)p s _(.- + ^)^i(^^ S ) 

[* [ { {a. - i/3)p s (a. - i/3)* \ 
+ J J M \Tr[(a.-i0)pU<r--il3)*] " P V W^K^^^-W)*]^^*)- 

Now, in order to address a diffusion approximation setup, we denote 

P = i\P\e l8 

and we investigate the limit e = v2/(\/7o|/3|) — ► 0. Physically, such a limit corresponds to a 
strongly excited local oscillator (|/?| — > oo). Naturally, here, we define 

Df = _ M a _ e -» + L and m= [l° a _ e -ie + i_ 
1 V 2 e 2 V 2 e 

It is then straightforward that 

_ 0,6-* _ (_ -. e -«)* + CT .e-^ - (a_e~ w y) =- (). 

This corresponds to the condition (|7|). Applying Theorem 1 and putting C = <7_e~ , we 
obtain the diffusion equation 

Pt = po+ / (-i[H,p s ]-^{a*c',pj-^{(7*c',p s } + ^c'p s c'* + ^ s c'*)d s 

J 

+ / -yffipP' + A'.C'* - Tr[ Ps (C + C*)]p s ) dWi(t) 

+ £ (Cp s + ^C* - Tr[p s (C + C*)]p.}dW 2 (t), (27) 

where ((Wi(i), W2(t)) is a 2-dimensional Brownian motion. Let us notice that we can define 
W t = yfl(W 2 (t)-W 1 (t)), which is a standard Brownian motion. The equation (I27p becomes 
then 

ft = Po + J L( Ps )ds + J ^To(Cp s + Ps C* - Tt[ Ps (C + C*)]p s yw t , (28) 
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where 

L(p) = -i[H,p] - ^{C*C,p s } + l0 Cp s C*. 

We recover the equation which has been presented, in the Introduction, as the model of 
homodyne detection. 

2.2 Diffusion Approximation in Heterodyne Setup 

Our last application is the model of heterodyne detection. In this case the parameter /3 is 
replaced by (3 *-* fit = f3e~ lI ^ t , where A is the detuning of the local oscillator. By assuming 
that the result of Theorem 1 is still valid for coefficients depending on time (there are no 
additional difficulties and the proof is the same). We get a similar expression for the het- 
erodyne detection. In this case the operator C is replaced by a time dependent operator 
C(t) = Ce lAt . Physically, a natural assumption is A > 1, that is A — > oo. Let us investigate 
how the equation (|28p is transformed under this condition. To this end, let us introduce some 
elements. We consider the following generators, defined for all / E C 2 , all p E S and all s E R 

by 

A A f(s,p) = D p f(L(p)) + ^D 2 p f(h(s,A,p),h(s,A,p)), 

Af(p) = D p f(L(p)) + ^D 2 p f(h + (p),h + (p)) + ^D 2 p f(ih^p),ih^p)), (29) 

where 

h(s,A,p) = ^(c(s)p + pC(sy-Tr[p(C(s) + C( S y)]p), (30) 
h±{ P ) = ^-(Cp±pC*-Tv[p(C±C*)]p). (31) 

Let po be a state and let (p^) and (pt) be the solutions of the problems of martingale associated 
with (*4a,Po) and (A,po). These solutions can be expressed as solutions of the following 
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stochastic differential equations 





o Jo 



Pt = 




o Jo Jo 



where (Wt), (W\(t)), and (W2(t)) are three independent Brownian motions, defined on a same 
probability space (£l,J-,P). Now, we are in the position to express the limit result. 

Proposition 3 Let (p^) be the solution of 13^) and let (pt) be the solution of $33\i . 

Therefore, the family of processes (pt) converges in distribution to the process (pt), when 
A goes to infinity. 

Proof: We will not show the tightness property of the family (p^) (one can show a similar 
result as Proposition 1 in the previous Section). Here, we cannot use of the convergence of 
generators as in Proposition 2. Actually, the convergence of Aa to A is not true since e jAi 
has no limit, when A goes to infinity. Here, we show directly that for all sequences A n , which 
converge to infinity such that (p t n ) converges in distribution, we have 



for all m > 0, for all < t% < t<i < . . . < t m < t < t + s, for all functions (#i)i=i,...,m and for 
all / in C7? (let us remind that in the previous section, this result was just a consequence of 
the convergence of generators). Let (A n ) be such a sequence and let (pt) be the limit process. 
Since (pf~") is valued is in the set of states, we first remark that the limit process (pt) is also 
valued in the set of states (the property satisfied by a state are closed for the topology of 



rt 



lim E 




(34) 
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Skorohod). Now, we can remark 



E 



f(pfrs)-f(P? n ) 



E 



t+s \ m 

' i=i 

t+S [t + S 

8=1 



t+S 



E 



i=l 



du. 



In order to estimate the limit of this term, we have to notice that 



(35) 



h(s,A,p) = V2(h + (p) cos(As) +ih-(p) sin(As)) . 



(36) 



for all p E S. It follows that 



A A J(s,p)-Af(p) = 

= -D 2 p f (y2(h+(p) cos(As) + ih^(p) sm(As)),V2~(h + (p) cos(As) + ih^(p) sia(As] 



\Dlf{h + {p)M{p)) + ^Dlf{ih-{p),ih^ 

l -D 2 p f{h + {p),h + (p)) (2cos 2 (A S ) - 1) + ~D 2 p f(ih„(p),ih-(p)) (2sin 2 (A S ) - l) 
+ 2D p (h + (p),ih-(p))cos(As)sm(As), 
{h + {p),h + (p))cos(2As) - ^L p 



1 

2 pJ 



1 D 2 J(h + (p), h + {p)) cos(2As) - \D 2 J(ih-(p),ih-{p)) cos(2As) 



+ D 2 (h + (p),ih^(p)) S m(2As), 



(37) 



for all p 6 S. Then, we get 



35D < 



t+S 



E 



rt+s 

+ L E 

[t+S 

+ l E 



-D 2 p f{h + {p^),h + {p^)))X{Oi{pt 



cos(2A n s)du 



1 \ m 

- J D 2 /(^_^),^_(^)) n^(pf- 



i=l 



cos(2A n s)<iu 



{Dlfih+ip^.ih-ipt)))^^) 



i=l 



sin(2A n s)cfoi. (38) 
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Let (i), (ii) and (Hi) denote the three terms of the sum. We have 



KOI < 



< 



E 



+ 



Re 



\Dlf{h + { P t)M{p% n )))tl e i(p£ n 

' i=l 

\D 2 p f{h + ^ u ),h + ^ u ))\\{e^ u ) 
i \ m 

~D 2 p f{h+{n u ),h+{v<u)))Y[Oi{n 



E 



du 



i=l 



e i2AnU du 



(39) 



Concerning the first term in this inequality, the term inside the integral f. s converges to 
zero, when n goes to infinity (according to the convergence in distribution of (p^") to (pt))- 
Moreover, this term is bounded independently of u and n since the processes (p^™) and (pt) 
take values in the set of states and the functions are C 2 . By applying the dominated Lebesgue 
Theorem, the first term converges to zero, when n goes to infinity. Concerning the second 
term, the term inside the integral in front of the term e lA " u is bounded and measurable (with 
respect to it), since A n goes to infinity, a classical result in Fourier transform theory, implies 
that the integral converges to zero. Let us stress that we can treat the other terms in a 
similar way and we show that the expression, defined by (|35|) . converges to zero when n goes 
to infinity and the result holds. 

Now, since the family of processes (p A ) is tight, we can apply a similar reasoning as the 
one used to prove the convergence in Theorem 1. Thus, for all sequences (A n ) converging to 
infinity, when n goes to infinity, we can extract a subsequence (A' n ) such that the property 
([M]) is satisfied for (p t n ). As a consequence the limit process of ((p t ")) is a solution for the 
martingale problem associated with (A, po). By a similar reasoning as the proof of Theorem 
1, we conclude that the family of processes (pt) converges in distribution to (pt) solution of 
the equation (j33"|) . □ 
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